Strong quantum interference in strongly disordered bosonic insulators 
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We study the variable-range hopping (VRH) of bosons in an array of sites witli short-range 
interactions and a large characteristic coordination number. The latter leads to strong quantum in- 
terference phenomena yet allows for their analytical study. We develop a functional renormalization 
group scheme that repeatedly eliminates high-energy sites properly renormalizing the tunnelling 
between the low-energy ones. Using this approach we determine the temperature and magnetic 
field dependence of the hopping conductivity and find a large positive magnetoresistance. With 
increasing magnetic field the behaviour of the conductivity crossovers from the Mott's law to an 
activational behaviour with the activation gap proportional to the magnetic field. 
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Transport in a huge variety of insulators is described 
by the model of variable-range hopping (VRH) [l[ based 
on a few rather general assumptions: charge carriers hop 
inelastically between random-energy sites, e.g. impuri- 
ties, grains, or (quasi-)localized states, the probability 
of a hop exponentially decreases with distance and has 
an activational temperature dependence. The particu- 
lar nature of the sites and the bath is not important, - 
the model can be applied to the transport of electrons in 
doped semiconductors as well as to that of Cooper pairs 
in macroscopically insulating superconductive materials. 

In the course of each inelastic hop between two sites 
a charge carrier may tunnel elastically via intermedi- 
ate sites. The interference between different-path contri- 
butions gives rise to mesoscopic and magneto-transport 
phenomena 0- The interference is suppressed in strongly 
insulating materials because there the elastic transmis- 
sion is dominated either by the direct tunnelling path 
between two sites or by the shortest chain of intermedi- 
ary sites. With increasing the intcrsitc couplings or the 
site concentration, more alternative paths come into play 
and the interference effects get more pronounced. 

This Letter is devoted to the study of the VRH trans- 
port in bosonic insulators with large coordination num- 
bers. In these materials the interference phenomena are 
particularly strong because of the constructive interfer- 
ence between nearly all tunnelling processes, including 
the many-body ones. Large coordination number Z ^ 1 
of the effective array of sites allows us to treat the in- 
terference effects analytically, unlike the previous VRH 
studies, and is characteristic of a broad class of insu- 
lating materials. Those include, for example, supercon- 
ductive films close to the metal-insulator transition (see, 
e.g., Ref. 3 for a review), densely packed granular ma- 
terials, disordered bosonic systems with many particles 
on the localization length, "superconductive hay" etc. 
Strong interference in these materials manifests itself, for 
instance, in a huge magnetoresistance magnetic field 
may change conductance by orders of magnitude or even 
drive a superconductor into a strong Cooper-pair insula- 



tor. Transport phenomena and superconductor-insulator 
transitions are also being observed since recently in artifi- 
cially made arrays of Josephson junctions with Z ^ iQ. 

In general, the interference in these materials is in- 
tractable perturbatively, unlike many earlier VRH mod- 
els (cf., e.g., Ref. 13). Tunnelling along paths of sites 
with close energies [Sf is also irrelevant for such systems 
due to the small portion of such paths among all possible 
ones. 

Model. Although we focus on the transport of bosons, 
the results of this Letter apply also to systems with an 
arbitrary charge carrier statistics if the VRH is domi- 
nated by single-particle hops. The latter may be the 
case, for instance, in an electronic system with a shallow 
Fermi sea, that is sufficiently exceeded by the charac- 
teristic Coulomb energy. As we consider systems with 
Z ^ 1, we also expect the results to hold qualitatively in 
the materials with coordination numbers of order unity. 

We assume for simplicity that the intersite interactions 
are short-range and the number of particles on each site 
uniquely determines their state. This simplification is 
possible if the other states are absent or separated by 
large energy gaps, which is characteristic of electrons 
in doped semiconductors l9||, Cooper-pairs in macroscop- 
ically insulating media[lClllll|. and of the charge carriers 
in nearly all other materials under consideration. 

The particle motion is governed by the Hamiltonian 

n=Y,Ur{nr)-Y.Ur'blbr'- (1) 
r r^r' 

Here C/r(?^r) is the energy of rir particles on site r, a large 
random quantity affected by disorder and stray charges 
close to the site. The last term in Eq. ^ accounts for 
the intersite hopping, 6j and br being (bosonic) creation 
and annihilation operators on grain r, rir = &J&r- 

Below we count the energies of all excitations from the 
energy of the "ground charging state" , the configuration 
of integers rir that minimize the first term in the right- 
hand side (rhs) of Eq. ^ . Typically, the energies and 
E~ , needed respectively to add or remove a particle from 



2 



the site r, are determined by the Coulomb interactions 
with the other particles and stray charges, and we refer to 
these energies and on-site excitation states as respectively 
"charging" energies and states. 

As only low-energy states are relevant for low- 
temperature conductance, one can conveniently rewrite 
the Hamiltonian in a reduced Hilbert space,- each site 
r has only 3 states: | ± l)r and |0)r, corresponding re- 
spectively to ±1 and extra particles as compared to 
the ground charging state. Introducing spin-1 operators; 
S^\ ± l)r = ±1 ± l)r, Sm = 0, S^\ T 1) = 7210), 
= V^l ± 1), we arrive at the reduced Hamiltonian 

^ = ^ E [i^r + E-) S^S^ + {E+ E-)S( 

r 

-^E^-'5+5-, (2) 

r.r' 

where the exact form of the coupling constants Jrr' de- 
pends on the microscopic details of the model. For in- 
stance, for a bosonic system with large occupancies rir ^ 
1 in the charging ground state, Jrr' = (nrnr')^/^trr'. 

In absence of inter-site interactions at strong disorder, 
the density of states (DoS) of the (anti)particle states has 
no Coulomb gapQ and, thus, reaches some finite value 
v at low energies. The assumption of a constant DoS in 
a certain energy interval is crucial for our consideration. 

Site-decimating RG. The transport of low-energy ex- 
citations can be analysed by constructing a logarithmic 
renormalization group (RG) technique that repeatedly 
eliminates the highest-energy charging states properly 
renormalizing the intersite tunnelling couplings between 
all the lower-energy states. 

Let us select a few (anti-)particle charging states with 
the very highest energies and split the Hamiltonian into 
the Hamiltonian T-Lhigh of the respective high-energy 
charging states, the Hamiltonian "Hiow of the rest of the 
system, and the tunnelling V between the selected states 
and all the others. Any eigenstate of the system is, 
strictly speaking, a superposition of the charging states 
on all sites. The projection '^low of the eigenfunction 
on the lower-energy charging states satisfies Schrddinger 
equation T-L'iow^iow = £'^iow with an effective (exact) 
Hamiltonian 

n'lon, = 'Hiow - VCHMgh - e)-^V, (3) 

£ being the eigenenergy. The latter may be neglected in 
Eq. dS]) if the energies of the selected states significantly 
exceed the energies of the other states. 

Assume, the intersite tunnelling elements are small, so 
that one can neglect (anti-)particle hopping between the 
highest-energy charging states. Then we find that the 
modified Hamiltonian ([3]) has the same form as the initial 



Hamiltonian ([2]) provided the couplings are renormalized, 

Jrr- ^ Jrr'+2^(rtRp+l) +2^ ^_ 

Rp Rh 

(4) 

where Rp and R/i label the eliminated particle and an- 
tiparticle states respectively. The increment of the mag- 
nitude of the coupling |Jrr'| is positive, as the energies 
and E^^ in Eq. ^ are positive by definition. 

Eq. (|4]) can be understood as follows. In addition to 
the direct tunnelling between sites r and r' a particle 
may cotunnel via virtual high-energy states in one of two 
ways: l)via an intermediate particle state on site R, or 
2)first create a particle-antiparticle dipole on sites r', R, 
then — on sites R, r. The respective high-energy particle 
or antiparticle state on the site R can be omitted pro- 
vided the coupling constants Jrr' arc renormalized. Since 
processes I) and 2) interfere constructively in a bosonic 
system, the renormalization increases the tunnelling. 

The generic Eq. (|4]) may be of little use for arbitrary 
distributions of the site energies and locations. However, 
in the case of a large coordination number Z the aver- 
aging over the site positions leads, as we show below, 
to a simple expression involving the averaged function 
Jrr'. Moreover, large coordination number suppresses 
the fluctuations of strongly renormalized couplings, since 
all the decimated states contribute with the same sign to 
each coupling, while the fluctuations result from multiple 
random-sign contributions. 

Repeatedly eliminating the highest-energy charging 
(anti-) particle states and introducing continuous coordi- 
nates in the coupling function we arrive at a functional 
RG equation with a running high-energy cutoff E: 

9/ Jrr' =XJ JrRjRr'rfR, I = \og{Ea / E) . (5) 

We assumed here that the quantity A — 2((nR) + l)v is 
independent of energy. 

On each RG step we considered the cotunnelling via 
only one intermediate site, which is justified for suffi- 
ciently small couplings J, 

AJ = y A Jrr' fir' < 1. (6) 

Indeed, taking into account the tunnelling between elim- 
inated charging states in the n-th order results in ~ 
[(AJ) log(£"/i?)]" relative correction to the tunnelling 
couplings on an RG step with energy rungs E' and E. 
Thus, Eq. ([5]) is the one-loop functional RG equation for 
the intersite tunnelling elements. 

Eq. ([5]) should be solved implying the coupling Jrr' 
corresponds to the bare one at the energies of order 
of the characteristic charging energy Eq in the non- 
renormalized Hamiltonian ([2]). The RG flow must be 
stopped if, as a result of the renormalization, the cou- 
plings become too large, violating condition or if the 
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high-energy cutoff E reaches either some characteristic 
energy of relevant conducting excitations or the mean 
level spacing on the characteristic radius ^ of the 

renormalized function Jj.^' in a d-dimensional space. 

The smallness of the fluctuations of the coupling on 
each RG step requires Z{E/ E')[\n{E' / E)]^ > 1. So long 
as the fluctuations are suppressed, the renormalized func- 
tion Jrr' fully describes the conduction. We emphasize 
that this is characteristic of a bosonic system, considered 
in the present Letter, and makes our analysis very differ- 
ent from that of fermionic insulators (il. [l^ . In a fermionic 
system our RG scheme would yield a non-renormalized 
coupling function with large fluctuations, which would 
determine then the transport properties of the system. 
The latter does apply, however, to a system with a so 
shallow Fermi sea that the transport is effectively single- 
particle,- the renormalization then occurs similarly to 
the bosonic case. 

Mean-field critical point. In absence of magnetic field 
the intersite tunnelling is translationally invariant, Jrr' = 
Jr_r'- Then for the Fourier-transorm Jk of the coupling 
function Eq. ([S]) yields 



Jk(0-^k(0)[l-AJk(0)r'- 



(7) 



where Jk (0) is the bare value of the tunnelling amplitude. 
In the limit of a strong insulator, AJk^ ^ 1, the renor- 
malization results in a small correction to the tunnelling 
and can also be obtained perturbativcly by considering 
the tunnelling processes via one intermediate site. 

The renormalization is more pronounced in less insu- 
lating materials. In particular, the renormalized tun- 
nelling elements diverge if the coordination number Z 
or the bare coupling are sufficiently large, namely, if the 
quantity Jk(0) exceeds some critical value Jc, 



XJJ r, 



= 1, 



(8) 



Imax being the maximal value of the logarithm that stops 
the RG flow. 

Eq. (lU matches the mean-field (MF) criterion of the 
phase transition in a system with Hamiltonian ^ and 
the uniform order parameter A = (S'r'") . However, we em- 
phasize that the above analysis is not sufficient to prove 
the existence of the metal-insulator or superconductor- 
insulator transition when condition ([8]) is fulfilled, be- 
cause we have to stop the RG flow if the inequality © 
no longer applies. Clearly, to study such transitions one 
has to employ alternative methods, which is beyond the 
scope of this paper. Now, we concentrate on the trans- 
port in the insulating regime when the effective ampli- 
tude Jk (0: Eq. dZ]), is stiU smaU. 

Conductivity away from the MF transition point. We 
can find rather easily the temperature dependency of con- 
ductivity in the region 



Eq. (j9]) means that the system is quite close to the MF 
transition and the effective tunnelling amplitude Jk (l) is 
strongly renormalized, although being small. Deeper in 
the insulating regime, i.e. at smaller J(0) than that in 
Eq. the renormalization is negligible and so is the 
effect of the cotunnelling on the VRH. 

Defining the characteristic radius C of the non- 
renormalized intersite coupling as 



Jk(0)« J(0)(l-C'k2), 







we find from Eq. ([7|) the dependence of the tunnelling on 
distance |r — r'| » (^i 

Jrr,(0- J(0)i. ' 



J(0) ^^-'^C^\r - r'p-'^expC-lr - r'|/0, (10) 



where we defined the characteristic radius 

1-1/2 



e = C[(AJ(0)/)-i-l] 



(11) 



of the renormalized coupling function. 

At sufficiently low temperatures the RG fiow should 
be stopped if the characteristic number Z{1) = Ev^'^ of 
states on length ^ approaches unity. Eq. (fTTj) and I = 
ln(£^o/^') define self-consistently the energy E = [vS,'^)^^ 
that terminates the flow. 

Thus, applying the RG procedure we arrive at a sparse 
array of sites with the coordination number Z of or- 
der unity, the characteristic intersite spacing Eq. (|lip . 
characteristic on-site energies Emin = 1/(!^C''), and the 
tunnelling amplitudes given by Eq. ([T0| . Inelastic hop- 
ping in the latter model is dominated by the direct hops 
between the sites, while the interference only slightly 
modifies the behaviour of the particle wavefunctions. 

Owing to the weakness of the tunnelling, each parti- 
cle excitation is almost entirely localized on one site. By 
construction of the site-decimating RG the particle am- 
plitude distance r away from this site estimates 

*(r) ~ Jo,MEa/Er)]/Er, Er = max[iiyr''y\Ern,n] 

(12) 

in the first order in the small tunnelling. At large 
distances the wavefunction decays exponentially cx 
exp(-r/,^), cf. Eqs. ^ and Applying the tradi- 

tional Mott's arguments [l| we find the temperature de- 
pendence of conductivity (Mott's law): 



ij = (To exp 



i/(d+iy 



T « 



KnL « 1 - J{^)/Jc « 1. 



(9) 



(13) 

with a renormalized Mott's temperature T'^j matching 
the value of the cutoff Emin m the RG procedure. 

At r > the RG should be stopped when the run- 
ning energy cutoff reaches temperature T, as the lat- 
ter sets the characteristic energy of conducting excita- 
tions. The power-lower hopping probability on distances 
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shorter than ^ leads to the temperature dependency of 
conductivity weaker than exponential. The standard 
Mott's arguments yield[l3] a power law a ~ T", although 
those cannot be used to find explicitly the power a. 

Effect of magnetic field. In magnetic field each tun- 
nelling amplitude acquires a phase: 



Jr 



Ur 



|cxp 



iqc 



A(r)(ir 



(14) 



Here q is the particle charge and the integration is car- 
ried out along a certain path connecting sites ri and 
r2, which in general depends on the microscopic details 
and for simplicity assumed below to be a straight line. 
We neglect the modification of the non-renormalized val- 
ues of the couplings |Jrir2l due to the "shrinkage of 
wavefunctions" [9| . In principle, it can be straightfor- 
wardly included in the magnetoconductance. 

The RG procedure leaves the phase factors intact but 
modifies the tunnelling amplitudes as 



A / I JrRlkR.r'1 CXp 



iqc 



A(r)dr 



Rr' 



dR. 
. 

The contour integral in Eq. (jisp gives the magnetic flux 
through the triangular contour r R — >■ r' — >■ r. 

Let us consider a two-dimensional system in a weak 
magnetic field B, so that the rcnormalization is not sup- 
pressed and the length 



Lb = [cl{qB)f^ 



(16) 



exceeds the characteristic coupling radius ^ in the ab- 
sence of magnetic field, Eq. pT|) . Then the RG should 
be stopped again at the same high-energy cutoff (i^^'^)^^. 

According to Eq. ([15]), the tunnelling on distance r 
is unaffected by magnetic field if the flux through the 
characteristic "interference area" r'^/^^-'^/^ is less than 
c/q, i.e. if r^l'^^l'^ < L\. Therefore, the magnetic field 
has no effect on conductance provided the mean inelas- 
tic hopping distance[l| r ~ ^(Ta^/T)^/^ does not exceed 
the scale R ^ L^'''^^^^/'^, which otherwise determines the 
characteristic length of the inelastic hopping, as it serves 
as a cutoff of the tunnelling Jrr' and the excitation wave- 
fimctions, Eq. ([T^ . 

At larger magnetic fields, ^ > Lb, the RG fiow should 
be stopped when the number of remaining states on 
length Lb approaches unity, and the renormalized cou- 
pling is cut on the length Lb-, Fig. [TJ The conduction 
comes then from the inelastic hopping on distance i? ^ 
Lb between sites of the characteristic energy l/(i/L^). 

Thus, we arrive at the conductivity 



a — exp 



1 min(e2/3,i^/') 



vT 



(17) 



that shows an activational behaviour with the activation 
gap oc _B in high magnetic fields and cx B'^/^ in lower 
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FIG. 1. The dependence of the renormalized tunnelling on 
distance, obtained from the numerical solution of Eq. (|15[) . 
Three curves correspond to different values of magnetic field 
[in units c/(gC2)] B= 1)0, 2)0.001, 3)0.01. I = 0.97^^^:.. 
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FIG. 2. (Colour online) The dependence of conductivity on 
magnetic field in a two-dimensional system. 



fields. Fig. [2] Upon decreasing the field so that T/T^j > 
(S./Lb)'^ the conductivity p7)) crossovers to Eq. ([T3| . 

At very strong magnetic fields, B ^ c/{qC'^), the 
interference between the tunnelling paths is destroyed, 
which suppresses the renormalization and leads to the 
usual Mott's-law conductivity, Eq. (fTB]). with = C- 

Recently a similar model of Cooper-pair hopping in an 
array of sites with random energies and large coordina- 
tion number Z ^ I has been studied in Refs. [3, [Tlj in 
order to describe the superconductor-insulator transition 
(SIT) in disordered superconductive films. In Rcf. [llj 
the effective array of sites has been mapped on a Bethe 
lattice and studied by means of an MF-typc approach 
with a site-dependent order parameter, using Z 3> 1. 
We disagree with such a mapping because Bethe lattice 
contains no finite loops of sites and, thus, cannot account 
properly for the interference phenomena considered here. 
Actually, the mapping of Ref. [ll| cannot be justified 
even after performing the rcnormalization, because the 
RG flow can be stopped only when the characteristic co- 
ordination number Z{1) reaches unity. 

An interesting issue that deserves a separate investi- 
gation is the role of the intersite interactions in the SIT 
in disordered superconductive films, as even arbitrarily 
small long-range charging interactions are known [l^ to 
rcnormalizc the parameters of a clean system and change 
the order of the phase transition. 

Conclusion. We considered low-temperature variable- 
range hopping transport of bosons in an array of sites 
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with random energies. We constructed a renormalization 
group technique that ahows one to repeatedly ehminate 
the highest-energy on-site states, renormahzing the inter- 
site couphngs between the lower-energy ones. Using this 
procedure, we reduced the system to a sparse array of 
sites and found the dependency of conductivity on tem- 
perature and magnetic field, Eqs. (fT3|) and p7)) . When 
increasing magnetic field the dependency crossovers from 
the Mott's law to Arrhenius law with the activation tem- 
perature proportional to the magnetic field. 
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